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Abstract 

Using effective formulas we analyze the Bose-Einstein correlations (BEC) data cor- 
rected for Coulomb interactions provided by STAR Collaboration and the quasi- 
corrected data (raw data with acceptance correction etc) on 27r and Svr BEC by us- 
ing Coulomb wave function with coherence parameter included. The corresponding 
magnitudes of the interaction regions turn out to be almost the same: i?coui(27r) ~ 
|i?Coui(3vr). i?coui means the size of interaction region obtained in terms of Coulomb 
wave function. This approximate relation is also confirmed by the core-halo model. 
Moreover, the genuine 3rd order term of BEC has also been investigated in this 
framework and its magnitude has been estimated both in the fully corrected data 
and in the quasi-corrected data. 



1 Introduction 



One of interesting subjects in heavy ion physics at high energies is the higher 
order Bose-Einstein correlation(BEC) as well as the normal (2nd order) BEC [1- 
3]. The investigation of the genuine 3rd order term of BEC is a current im- 
portant topic. We are also interested in this subject to examine usefulness 
of a theoretical formulas described by Coulomb wave functions. Moreover, it 
is well known that information on the magnitude of the interaction region is 
necessary to estimate the energy density of produced hadrons [4]. Actually 
several Collaborations have reported interesting results in Refs. [5], [6] and 
[7]. 
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Very recently the STAR Collaboration [8] has reported new data on BEC 
for three negatively charged pions (Svr^), with Coulomb corrections imposed, 
where i? = 5 fm is assumed as input [5,8-10]. Their data on 2tt BEC have 
been already published in Ref. [10], whereas the quasi-corrected data (raw data 
with acceptance correction etc) on Stt" BEC have been shown in Ref. [9]. In 
Table 1 we provide our classification of different kinds of data available [11- 
14] ^ . We shall be interested here in different kinds of data provided by the 
same experiment, to know description of them. We shall be able to examine 
whether our formulation of BEC in terms of Coulomb wave function as given 
in Ref. [14] is useful when apphed to concrete data. 

Table 1 

Category of data. 

1) Raw data raw data 

2) Quasi-corrected data (Q-CD) (raw data) x i^acceptance 

3) Gamow corrected data (raw data) x -K^acceptance 

X Kg amow 

4) Coulomb corrected data (raw data) x ^acceptance x -f^Couiomb 
~Z . , , , , ~ _ JUd'^Xip{xi) ■ [plane w. f.|^ 

-'^Coulomb IS CaiCUiateu as -fvCoulomb — rr-Tj'? r \ i i ri9' 

J Y[d'^Xip{xi) ■ |Coulomb w. i.p 
where p{xi) = {fl'^li^fl'^ &yi-^{-0^xf[ [11-13]. Notice that R = input 

(/? = l/-\/2i?) is necessary for -f^Couiomb) where i? = 5 fm is used in 

Refs. [5,8-10]. For the systematic error bars we assume 5% for Sir BEC 

and 0.5% for 27r BEC. 



For various kinds of data in Table 1, we prepare following two methods^ . 

i) A conventional method: Effective formulas derived from the plane wave 
with the source function, p{xi) are applied to 4) Coulomb corrected data 
categorized in Table 1. 

ii) A method proposed in Ref. [14]: A formula described by the integration 
of the Coulomb wave function with p{xi) is applied to 2) quasi-corrected 
data(Q-CD) in Table 1. This will be explained in the next paragraph. 

Let us start with STAR data including Coulomb corrections as given by [8,10]. 
We analyze these data using the following effective formulas, introducing pa- 
rameters, i?eff, A2 and A3, 



C2 = l + A2e-(^«««i2)% (1) 

^ The aTithor of Ref. [11] pointed out that the Coulomb correction is more accurate 
than the Gamow correction. 

^ Different methods for data of 27r BEC in Refs. [15,16] will be explained in the 
third paragraph. 
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(2) 



where Q^. ^ (k. _ ^.)2 ^nd = (ki - + (ks - kg)^ + (kg - ki)^. 

Our estimated parameters for 27r~ and Svr" BEC are shown in Table 2. To 
reproduce R ~ 6.3 fm reported by STAR Collaboration, we have to skip first 
two data points, see Fig. 1(b). 
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Fig. 1. (a) Analysis of data [10] by Eq. (1). (b) Data with (0) are skipped in analysis, 
(c) Analysis of data [8] by Eq. (2). 



Table 2 

Data (Coulomb correction with ii = 5 fm (input)) on 27r and Stt by STAR Collab- 
oration [8,10] and our analyses of data by Eqs. (1) and (2). 



Data 


Theor. 


C 


RcS [fm] 


A2 or A3 




27r [10] 


Eq. (1) 


0.999±0.003 


7.84±0.20 


0.62±0.02 


30.8/24 


27r [10] 


Eq. (1) 


0.996±0.003 


6.52±0.31 


0.43±0.03 


8.32/22 


Stt [8,9] 


Eq. (2) 


0.999±0.012 


7.83±0.26 


0.57±0.02 


5.32/35 



On the other hand, in Ref. [14] we have proposed another theoretical formu- 
lation of the Coulomb wave function approach to BEC containing the degree 
of coherence parameter A. We shall explain it in the next paragraph in de- 
tail. Here it will be applied to the raw data (with the acceptance corrections 
etc included) and its result will be then compared with the previous one in 
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paragraph 3. Therein we shall analyze also the raw data on 27r~ BEC. In 4th 
paragraph, we are going to analyze the data by the core-halo model [13,14]. 
Our concluding remarks will be in the final paragraph. It will also contain dis- 
cussion of the problem of phases among the three 2nd order BEC correlation 
functions C2, which will be done by introducing a parameter A (which can 
attain complex values). Finally we compare our results with those obtained 
by means of formula given in Ref. [17]. 



2 A theoretical formula for charged Stt BEC. 



We shall provide now brief explanations of our formulas (sec Ref. [14] for more 
details). The charged 3n~ BEC should be described by the following products 
of three Coulomb wave functions, see Fig. 2. 
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Fig. 2. Diagram reflecting three-charged particles Bose-Einstein Correlation (BEC) 
and Coulomb potential (K). X means the exchange eflFect of BEC. 



The theoretical formula for 37r BEC is expressed as 



4/n 

" •' i=l 



V'kika (Xl, X2)V'^,k3(X2, X3)V'£ki (^3, Xi) 

+^^ik2 (Xl, X3)z^^^k3(x3, X2)V^k3ki (X2, Xi) 
+^kik2(x2, xO^/'gkslxi, X3)^/'k3ki (X3, X2) 
+^kik2(x2, X3)^/'£k3(x3, Xi)^g^ki(xi, X2) 
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+V'kik2(x3, Xi)V'k2k3(Xl, X2)V'k3ki(x2, X3) 



C 



+'0kik2(x3,X2)V'k2k3(x2,Xi)^^3k,(xi,X3) (3) 



c 



Here ^kikj (xi, x^ ) are the Coulomb wave functions of the respective 2-body 
coUision expressed as, 



^gk,(xi, X,) = r(l + ir;,,-)e"^-/'e^'^----F[-ir7,,-, 1; i(%r,, - k,,- • r,,)], (4) 

with Yij = (xj-Xj), kij = {ki-kj)/2, = \rij\, kij = \kij\ and r]ij = ma/kij. 
F[a, b; x] and r(a;) are the confluent hypergeometric function and the Gamma 
function, respectively. 



In order to introduce a parameter of the degree of coherence (A) in Eq. (3), 
we decompose the Coulomb wave functions assigning to amplitudes A{i) the 
following Coulomb wave functions (taken in our calculation in the plane wave 
approximation (PWA)): 



^(l)=^kik2(Xl' X2)^/'k,k3(x2, X3)t/'g^ki(x3, Xi) 

PWA giki2-ri2gik23-r23gik3i r3i _ g(3/2)j(ki ■xi+k2-X2+k3-X3) 

^(2) =^^ik2(xi, X3)?/'^2k3(x3, X2)^£ki(x2, Xi) 

gjki2-ri3gjk23 r32gik3i-r2i _ g(3/2)i(ki •xi+k2-X3+k3-X2) 

^(3)=V'kik2(x2, Xi)V'k2k3(Xl, X3)V'k3ki(x3, X2) 

gjki2-r2igik23Ti3gik3i-r32 _ g(3/2)i(ki-X2+k2-xi+k3-X3) 

^(4)='0kik2(x2, X3)V'k2k3(x3, Xi)V'k3ki(Xl, X2) 

PW^ gjki2T23gjk23T3igik3i-ri2 _ g(3/2)i(ki-X2+k2-X3+k3-xi) 

^(5)='0kik2(x3, Xi)V'k2k3(Xl, X2)V'k3ki(x2, X3) 

PW^ giki2-r3igjk23-ri2gik3i-r23 _ g(3/2)i(ki •X3+k2-xi+k3-X2) 

^(6)=<k2(x3, X2Xk3(x2, Xi)^g^ki(Xl, X3) 

PW^ gjki2-r32gik23-r2igik3i-ri3 _ g(3/2)i(ki •X3+k2-X2+k3-xi) 

Using A{i),^ we obtain the following three groups denoted as Fj's, 



^ The normal plane wave amplitude ^piane(^) are expressed as 

A , n — P*('^l'''l+'^2-X2+k3-X3) A /'9\ _ J(ki-xi+k2-X3+k3-X2) 

See Ref. [2]. 



(5a) 
(5b) 
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F,^lY.mA*it)'^l, (6a) 

F2 = -[A{1)A*{2) + A{1)A*{3) + A{1)A*{6) + A{2)A*{4) + A(2)A*{5) 
6 

+A{3)A*{A) + A{3)A*{5) + A{A)A*{6) + A{5)A*{6) + c. c] (6b) 
BEC between two-charged particles (6c) 

F3 = ^[A(1)A*(4) + A(1)A*(5) + A(2)A*(3) + A{2)A*{6) + A{3)A*{6) 

+^(4)^*(5) + c. c] (6d) 
BEC among three-charged particles (6e) 

We have to assign the degree of coherence ( real number of \/A) to the cross 
marks (x) in Fig. 2. We obtain finally that 



Cjfl d'x^pi^^) ■ [Fi + AF2 + X'^'Fs] , (7) 



where p(x) = {p^/nf/^ ex.p[-p'^x^] {(3 = l/V^R) and A is the degree of 
coherence. Moreover, we can also derive simple formula which uses only plane 
wave approximation (PWA). It can be written as 



1^ 



^ Eq. (7) ^laV I 1 + 3A E e-?«^«^. + 2X'/'e-i^'^'^ \ , (8) 

>j 



Actually, Eq. (8) is also applied for the Gamow corrected data, because it 
contains no Coulomb effect. R stands for i^camow 



3 Analyses of data on BEC. 



a) Analyses of data on Stt" BEC by Eqs. (7) and (8): In Ref. [9] quasi- 
corrected data at y^s/vw" = 130 CeV have been reported. We can confirm them 
by using method proposed in Refs. [13,14] and results are shown in Table 3 
and Fig. 3. Eq. (7) and Eq. (8) with the Gamow correction i^camow are used 
in the present analyses. It is very interesting that the corresponding values of 
size parameter R estimated by Eqs. (7) and (8) are almost the same. 

b) Analyses of data on 27r^ BEC: In Ref. [10] quasi-corrected data (Q- 
CD) on 27r~ BEC are given, which we have analyzed using following theoretical 
formula: 
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Fig. 3. (a) Analysis of quasi-corrected data (Q-CD) on Stt^BEC [8] by Eq. (7). (b) 
Analysis of quasi-corrected data (Q-CD) on 27r~ [10] by Eq. (9). 

Table 3 

a) Analyses of two kinds of data on Stt" BEC. R means i?Coui or ^Gamow b) 
Analyses of two kinds of data on 27r~ BEC by means of Eqs. (9) and (10). 



a) 3 TT" BEC Data 




R [fm] 


A 


xVn-d.f. 


quasi-corrected data(Q-CD) 

(Q-CD)/i^Gamow 


Eq. (7) 
Eq. (8) 


5.34±0.24 
5.03zb0.20 


0.56±0.02 
0.61±0.03 


2.80/35 
6.41/35 


b) 2 TT- BEC Data 




R [fm] 


A2 


xVn.d.f. 


quasi-corrected data(Q-CD) 

(Q-CD)/irGamow 


Eq. (9) 
Eq. (10) 


8.75±0.31 
7.41±0.20 


0.58±0.02 
0.61±0.02 


23.0/25 
39.4/25 



jY(2+ or 2-) 



^BG 



C j nd=^x,p(x,)-[G'i + AG'2], 



(9) 



where 



and 



V'kik2(Xl, X2) '+ '0kik2(x2, Xi) 



G'2 = Re«*k^(xi, X2Xk,(x2, xi)) . 
Applying the plane wave approximation (PWA) to the above equation ^ , we 

^ It is worth while to mention a different method for 27r BEC utilized in Refs. [15] 
based on theoretical studies of Ref. [16]. The following formula is applied to 2) 
Quasi-corrected data (Q-CD) including the momentum resolution in Pb -|- Au col- 
lision at 40, 80 and 158 AGeV [15], 

C2 = C[1 + A'((1 + G)F*-1)]. 

Here G = eyiY)[—R\q\ — RoQo ~ ^'s^'s ~ '^^hilL^lo]- Rl, Ro, Rs and cross-term 
RoL are the Gaussian source radii, whose source functions are assumed, qf^^ = 
(ki - k2)2 - {El - E2f, and F* = w{kt) ■ (Fcoul(gmv) - 1) + 1, where Fcoul(ginv) 
denotes the Gaussian integration of the Coulomb wave function times the source 
function with R = 5 fm. See Ref. [15] for the notations, where the Monte Carlo 
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obtain the following expression, 



C2/i^Gamow = [1 + A^g-^^^^^-)'] , (10) 

where i?Gamow = RcS (cf., Eq. (1)). The results are also given in Table 3. It is 
worth to notice that estimated values of the size parameter R coincide with 
each other. Actually, in case of plane wave approximation (PWA) they are 
reduced to the following value: 



3 

2-Rcoul(37r) ~ -Rplane(37r) fti -Rplane(27r) (11a) 

3 

= - X 5 fm 7.5 fm (lib) 

The interaction region of i?piane(27r) and it!piane(37r) may be explained by the 
relation: 

R{Au) = 1.2A^ (Au) ~ 7 fm . (12) 



4 Analysis of data by core-halo model 



To confirm the results obtained previous paragraph, we use the core-halo 
model proposed in Ref. [2,13,14]. Introducing parameters fc = (^core)/('^tot) 
(the fraction of multiplicity from the core part), Pc — ('t-co) / (^core) (the fraction 
of coherently produced multiplicity from the core part), and p — (richao) / (^core) 
(chaoticity parameter defined in the laser optical approach), we have the fol- 
lowing expressions with p = 1 — pc, 



Ar(3-) 



Ar(2-) 



C 



C 



n d^x,p,{^,) ■ (Fi + /,VF2 + /.VFs) 

i=l 

+ / n c?'Vc(x.) • 5\-K2) (2/2p(l - p)F2 + 3/,V(l - p)F, 

i=l 

^ 1=1 



, (13) 



method is also introduced. In this approach the effective degree of coherence A = 

X'w{kt) • (i^Coul(ftnv) — 1) is expressed by the parameter A', the weight function 
w{kt) and Fcoulilmv), where kt = ^\^it + ^2t\- Finally, physical information on Rl, 
Rq, Rs and RoL are obtained as function of kf. Various theoretical formulas with 
w{kt) = 1 are investigated in Ref. [16]. 
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+ 



/ lld'xM^,) ■ 53(x2)2/2p(l -p)G2 
1=1 



(14) 



Here pc stands for the source function of the core part. The halo part is 
reflected by the delta function. 

Our results are shown in Table 4 and Fig. 4. As seen in them, the interaction 
ranges of R (core part) are estimated in the ranges of 5.3 fm < -Rcoui(37r) < 
7fm and 8.7 fm < Rcouii'^'^) < Hfm, respectively. The common region be- 
tween BEG of 37r~ and 27r~ is roughly described by 0.6^p^l.O and 0.75^/c^0.85. 
Actually we have confirmed Eq. (11), 1.5i?coui(37r) ~ i?coui(27r). 

Table 4 

Typical results from analysis of data by STAR Collaboration in terms of Eqs. (13) 
and (14). The effective degrees of coherence are defined by A3 = f^ip^ + 2p(l —p)) + 
fciP^ + 3p^(l — P)) and A2 = fcip^ + 2p(l — p)), respectively. Note that figures at 
p = 1 reproduce previous results in Table 3. 



P 



1.0 



0.8 



0.6 



fc 
R (fm) 

fc 

R (fm) 



Stt BEC 
0.75±0.02 0.78±0.02 



0.97 
5.34±0.24 

2.80/34 



1.00 
5.99±0.28 

1.39/34 



2tt BEC 
0.76±0.01 0.78±0.02 

0.58 0.59 
8.75±0.31 10.16±0.38 
23.0/24 16.3/24 



0.85±0.02 

1.00 
6.48±0.32 

0.98/34 

0.83±0.02 

0.58 
11.12±0.42 
16.1/24 



5 Concluding remcirks 



Because quasi-corrected data at ^/Jnn = 130 GeV have been reported by 
STAR Collaboration, we have applied theoretical formula, Eq. (7) to analyze 
them. Our method and the conventional method for the charged BEC are 
explained in paragraph 1. From our results, it can be said that both methods 
are (almost) equivalent to each other (when accounting for the Fourier factor 
(3/2) in Eq. (5) in case of charged 37r^ BEC). From our present study we 
have known that magnitude of the interaction region of Au + Au collisions 
at ^/sNN — 130 GeV is described by the relation it!coui(27r) ~ |-Rcoui(37r) 
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0.6' ' ^ ' ^ ' ^ ' ^ ' 1 

0.2 0.4 0.6 0.8 1 
P 

Fig. 4. Sets of fc and p estimated in analyses of BEC of 27r~ and 3tt~ . Widths of 
/c's stand for error bar of ±2(7. Common region are marked by mesh of a net. 

whose magnitude is explained by -R(Au) = 1.2 x A^^^ fm. This fact is also 
confirmed at p = 1 in Table 4, provided that the core-halo model proposed in 
Refs. [2,13,14] is used. This kind of study should be applied also to other data 
at SPS and RHIC energies [18]. 

Finally we consider problem of the phase factors appearing between three 2nd 
BEC C2(ki, k2), C2(k2, ka), and C2(k3, ki). According to Rcf. [17], we analyze 
the respective data by the following formula for the 3rd order BEC, 



Cs = C[l + X3YI e-(^«««»^)' + 2A^/'e-°-'(^««'5^)' x W] , (15) 

i>j 

where W = cos(0i2 + 023 + 4>3i), which, in the simplest form, is parameterized 
as W = cos{g x Q3). 

Similar formula can be obtained by replacing in Fig. 2 and Eq. (8) degree of 
coherence parameter \/A by three complex numbers (\/Ae*'^): 



yY(37r ) „ 3 

= C y n d^'x.pi^i) ■ [F, + AF2 + A^/^Fs X W] . (16) 

Our results are shown in Table 5 and Fig. 5. The genuine 3rd term of BEC 
in the quasi-corrected data (Q-CD) is explicitly shown in Fig. 5. It should be 
stressed that the phase factor is playing important role in formula N^^'^ /N^^ [19] , 
leading to smaller minimum of x^i see Table 5. The smaller x^'s in Table 5 
than those of Table 4 suggest us the usefulness of Eq. (16). To elucidate this 
fact, study on g is also necessary in a future [18] ^ . 

^ There is a possible interpretation as 5 = 30 GeV~^ ~ 6 fm, provided that 0.2 GeV- 
fm = 1 is used. 
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Table 5 

Analyses of 3it~ BEC by means of Eq. (15) and Eq. (16) with W = cos{g x Q3). 



Data 




R [fm] 


A3 or A 


9 [GcV -1] 


xVn.d.f. 


Coulomb CD 


Eq. (15) 


7.08±0.42 


0.62±0.03 


32.9±6.2 


1.00/34 


(Q-CD) 


Eq. (16) 


4.91±0.36 


0.61±0.04 


31.9±8.8 


0.57/34 




-0.05 I 

50 100 150 200 



Q3 (MeV) 

Fig. 5. (a) Analysis of Q-CD by Eq. (16). (b) The magnitude of the genuine term, 
(c) Behavior of cos{g x Q^). 
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